%This script computes the geomagnetic field values from Pomme
% plots the contours on Indian region

%20-6-2005

fday = tmjd(2005,1,1,0,0,0,1);
tha = [8:15];
phi = [74:81];
alt = 0; % at msl
B = get_pomme(fday,tha,phi,alt);
K = squeeze(B(3,:,:));
[x,y] = meshgrid(phi,90-tha);
%[map,ml] = gtopo30('c:\manoj\open\topo\E060N40',10,[4,17],[72,83]);
%worldmap(map,ml,'dem');
%s = worldhi([0,20],[70,90]);
%linem(s.lat,s.long);
worldmap('World');
%[c,ch] = contourm(flipud(y),x,K,[20000 15000 10000 5000 0 -5000 -10000],'color','k');
[c,ch] = contourm(flipud(y),x,K);
clabelm(c,ch);
%clabelm(c,ch,'manual');
set(ch,'FaceColor','none');
set(ch,'EdgeColor','k');
%gridm;
hidem(gca);
%[c,ch] = contourfm(flipud(y),x,abs(K),[5000 0]);
